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I. Introduction 


This manual contains descriptions of EASY5 macro component 
models developed for the spacecraft power system simulation 
sponsored by NASA GSFC. A brief explanation about how to use 
the macro components with the EASY5 Standard Components to 
build a specific system is given through an example in Sec- 
tion III. 

The macro components are ordered according to the follow- 
ing functional group. 

1. Converter power stage models 

2. Compensator models 

3. Current- feedback models 
. 4. PWM models 

5. Load models 

6. Solar array models 

7. Shunt regulator models 

The format of the component model descriptions is similar 
to the format of the EASY5 Standard Components descriptions 
in the EASY5 User's Guide. The circuit model of each macro 
component model is included along with some key equations and 
the model program is also attached for each model. The port 
inputs and outputs used for interconnection among models are 
shown on the right and left side of the box with arrows. 
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II. Macro component descriptions 


This section describes a set of EASY5 macro components 
developed in this study. Major equations, a circuit model 
and a program listing are provided for each macro component. 
The macro components are listed in the following table. 


TABLE OF MACRO COMPONENTS 


Macro 

Component 

Name 

Description . • 

Page 

BC 

- CONVERTER POWER STAGE MODELS - 
Buck Converter Power Stage 

4 

BT 

Boost Converter Power Stage 

8 

FB 

Buck/Boost Converter Power Stage 

12 

FW 

■ Foward Converter Power Stage 

16 

BP 

Battery-discharger Power Stage 

20 

DG 

- COMPENSATOR MODELS - 
Lead- lag Compensator 

22 

MP 

Two-pole One-zero Compensator 

24 

PZ 

Two-pole Two-zero Compensator 

26 

SM 

- CURRENT -FEEDBACK LOOP - 
SCM Current Loop 

28 

CC 

CIC Current Loop 

30 
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Macro 

Component 

Name 


Description 


Page 


- PWM MODELS - 


WM 

PWM (Constant Frequency Control) 

32 

FF 

PWM (Constant Off -Time Control) 

34 

NN 

PWM (Constant On-Time Control) 

36 


- LOAD MODELS - 


LO 

LCR Load 

38 

PT 

Constant Power Load 

40 


- SOLAR ARRAY MODELS - 


AR 

Solar Array 

42 

AS 

Solar Array Switching Unit 

46 


- SHUNT REGULATOR MODELS - 


HA 

Type-1 Shunt Regulator (Transfer- Function Model) 

52 

HB 

Type-1 Shunt Regulator (State-Equation Model) 

57 

FS 

Type-2 Shunt Regulator 

61 

AP 

Solar Array/Partial Shunt 

63 






BUCK CONVERTER POWER STAGE 


BC 



VI 

variable 

input voltage 

Volts 

12 

variable 

input current from load 

Amps 

IQ 

variable 

switching function 




(if switch is on, IQ=1) 




(if switch is' off, IQ=0) 


c 

parameter 

capacitance 

Farads 

L 

parameter 

inductance 

Henries 

RC 

parameter 

capacitor effective resistance 

Ohms 

RL 

parameter 

inductor effective resistance 

Ohms 


OUTPUT 


V2 

1 7 ^ T 

variable 

output voltage 

Volts 

11 

variable 

output current to source 

Amps 

IL 

state variable 

inductor current 

Amps 

VC 

state variable 

capacitor voltage 

Volts 

IS 

variable 

switch current 

Amps 

VL 

variable 

inductor voltage 

Volts 






EQUATIONS 



For IQ = 0, IL > 0 : 


♦ 



For IQ = 0, IL = 0 : 





ft*##***#**#####*#**#**#*#***#**#** 

** LARGE SIGNAL MODEL OF BUCK(W/Q LOAD) 

* POWER STAGE ONLY 

•ft*******#**#***#*****#*##*#*##***#****#* 

** REVISED 12/3/85 
MACRO FILE NAME"MACROS 
DEFINE MACRO = -BC 
MACRO INPUTS^ 

C L. RL RC 

IQ VI 12 

* IQ 5 SWITCHING FUNCTION 

* IQ=1 (SWITCH; ON) 

* IQ=0 (SWITCH; OFF) 

* VI ; INPUT VOLTAGE 

* 12 ; INPUT CURRENT (FROM LOAD) 

MACRO OUTPUTS- 

IS IL 

ILL VL 

VC V2 II 

* IL ; INDUCTOR CURRENT (STATE) 

* VC ; CAPACITOR VOLTAGE (STATE) 

* IS ; SWITCH CURRENT 

* ILL; INDUCTOR CURRENT (DUMMY) 

* VL ; TRANSFORMER VOLTAGE 

* V2 ; OUTPUT VOLTAGE 

* II ; OUTPUT CURRENT (TO SOURCE) 

MACRO CODE 

MACRO STOP SORT 

IF ( DABS( IQ BC -— ) . LT. 1 . E-10)THEN 
IF( IL BC — . LE. 0. ) THEN 
ILLBC — * O. 

MACRO DERIVATIVE, IL BC— =0. 

GOTO +++ 77 
ELSE 

**** TOFF ***** 

A1 1=- (RL BC — +RC BC — ) /L BC — 
A12=-l/L BC- ■ 

A21 = l/C BC — 

A22— O. 

B 1 1 =0 

B12=RC BC— /L BC — 

B21-0. 

B22=-l /C BC — 

C11=RC BC — 

Cl 2=1 
C21=0 
C22=0 
Dll =0 

D12=-RC BC — 

D2 1 =0 
D22=0 

* 

IS BC— =0 
END IF 
ELSE 

#** TON *** 

A11=-(RL BC — +RC BC — )/l. BC — 
A12=-l/L BC — 

A21 = l/C BC — 
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A?2~0 

Bll = t/t. BC — 

B12=RC IK: - -/L BC — 
B2 J =0 

B22=-l/C BC— 
CU=RC BC — 

Cl 2=1 
C21 = l 
C22=0 
D1 1=0 

D12=-RC BC — 

D21=0 

D22=0 

IS BC — =IL BC— 

END IF 
** OUTPUT EO. 

* 


- - 


*! V2 


Cll 

C12 

IL 


) DU 

D 12 ; 

5 VI ! 

*i 

*i 11 
«• 


C21 

1 

C22 1 
1 

VC 

+ 

1 D21 
1 
1 

1 

D22 ! 
1 
I 

I 12 j 


* 

V2 BC— =C11 *IL BC — +C12*VC BC — 

& + D11*V1 BC— +D12*I2 BC- 

Il BC — =C21*IL BC — 

VL BC — =<A11*IL BC— +A12*VC BC H311*V1 BC— + 

& B 12*12 BC— >*L BC — 

** STATE EG. 


* 


if 

if 

dlL/dt 

! All 

A 12 i 

t 

IL 

:bu 

+ i 

B 12 

5 VI 
1 
> 

if 

if 

: dVC/dt 

i 

! API 
« 

A22 : 

f 

i 

VC 

IB21 

1 

B22 

! 1.2 
t 


* 

MACRO DERIVA‘1 1VE. IL BC — =VL BC — /L BC — 

+++77 CONTINUE 

MACRO DERIVATIVE. VC BC — =A21*IL BC t-AP2*VC BC — * 

& B21*V1 BC — +B22*I2 BC — 

* 

ILLBC- — ~IL BC— 

IF < ILLBC— . LE. 0. ) ILLBC — =0. 

MACRO RESUME SORT 
END OF MACRO 
MODEL DESCRIPTION 
LOCATION= 22. BC 
END OF MODEL 
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BOOST CONVERTER POWER STAGE 


BT 



O-JvW^YV— ■ 
/ 

H X- r 0 

h 

x 2 » 

1 -j 

L 

► 

V L 

IQ 

c *1 

L. 



► 


INPUT 


VI 

variable 

f 

input voltage 

Volts 

12 

variable 

input current from load 

Amps 

IQ 

variable 

switching function 
(if switch is on, IQ=1) 
(if switch is off, IQ=0) 


c 

parameter 

capacitance 

Farads 

L 

parameter 

inductance 

Henries 

RC 

parameter 

capacitor effective resistance 

Ohms 

RL 

parameter 

inductor effective resistance 

Ohms 


OUTPUT 


V2 

variable 

output voltage 

Volts 

11 

variable 

output current to source 

Amps 

IL 

state variable 

inductor current 

Amps 

VC 

state variable 

capacitor voltage 

Volts 

IS 

variable 

switch current 

Amps 

VL 

variable 

inductor voltage 

Volts 
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EQUATIONS 
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ft** ft ftftftftft ft *«**•#*« * * ft ft ftftftftftft ftftftft ft ft * 

ft* LARGE SIGNAL. MODEL OF BOOST (W/O LOAD) 

* POWER STAGE ONLY 
ft REVISED 11/1 9/85 

ft^ftftftftftftftftftftftftftftftftftftftftftftftftftftftftftftftftftftftftftftftftft 

MACRO FILE NAME-MACROS 
DEFINE MACRO- BT 
MACRO INPUTS- 

C L RL RC 

VI 12 

MACRO OUTPUTS* 

IS IL 

ILL 

VC V2 11 

MACRO CODE 
MACRO STOP SORT 

IF (DABS( IQ BT— ).LT. 1. E-5)THEN 
IF < IL BT — . LE. O. )THEN 
ILLBT — -O. 

MACRO DERIVATIVE, IL BT — =0. 

GOTO *++77 
ELSE 

ftftftft TOFF ftftftftft 

A11=-<RL BT — +RC BT — )/L BT— 
A12=-l/L BT — 

A21=l/C BT— 

A22=0. 

B 1 1 = 1 /L BT — 

B12=RC BT — /L BT — 

B21=0. 

B22=-l/C BT— 

C 1 1-RC BT- - 
Cl 2=1 
C21 = l 
C22=0 
D1 1=0 

D12=-RC BT — 

D21=0 

D22=0 

ft 

IS BT— = 0 
END IF 
ELSE 

*** TON *** 

A1 1=-<RL BT — )/L BT — 

A 12=0 
A2i=0 
A22=0 

Bll=l/L BT — 

B 12=0 
B21=0 

BS3——1 / C BT — 

Cl 1=0 
Cl 2=1 
C21=l 
C22=0 
Dll =0 

D12=-RC BT — 

D21=0 
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D22-0 

IS BT — = n B I — 
END IF 
Q BT— = 10 
** OUTPUT EG. 

«• 

* — . _ 


V2 


Cll 

Cl 2 ! 
1 

! IL 
» 


i Dll 
1 

D12 i 
1 

| VI 

11 


C21 

CP2 ! 

s vc 

T 

! D21 

D22 ! 

! 12 


x. II II I I II 

* I I II I I II 

* 

V2 BT — =C11 »IL BT — +C12*VC BT — 

& + D11*V1 BT— +D12*I2 BT— 

II BT — -C21*IL BT — 

VL=(A11*IL BT — +A12*VC BT — +B11#V1 BT — + 
& B12*I2 BT — >*L BT- - 

** STATE EG. 

* 


* 

dlL/dt i 

All 

A 12 i 

i IL 

IB .11 

b 12 : 

i VI ! 

K 

I 

> 

1 

i 

+ i 

* 

1 

I 1 

1 » 

tt 

dVC/dt | 

! A 21 

AP 2 ! 

i vc 

! B 21 

B 22 : 

i 12 ; 

■* 

* 

1 

1 

1 

j 

i 

f 

1 

1 

1 

i » 

» > 


* 


MACRO DERIVATIVE, IL BT— «VL/L BT — 

<-++77 CONTINUE 

MACRO DERIVATIVE, VC BT— =A21*IL BT — +A22*VC BT — + 
& B21*V1 BT— +B22*I2 ST- 

S’ 

ILLBT --" IL. BT — 

IF ( ILLBT — . L.E. 0. JILLBl — -O. 

Vf 

MACRO RESUME SORT 
END OF MACRO 
MODEL DESCRIPTION 
LOCATION- 22, BT 
END OF MODEL 
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FB 


BUCK/BOOST (FLYBACK) CONVERTER POWER STAGE 



INPUT 


VI 

variable 

input voltage 

Volts 

12 

variable 

input current from load 

Amps 

IQ 

variable 

switching function 
(if switch is on, IQ=1 ) 
(if switch is off, IQ=0) 


c 

parameter 

capacitance 

Farads 

LP 

parameter 

primary inductance 

Henries 

LS 

parameter 

secondary inductance 

Henries 

RC 

parameter 

capacitor effective resistance 

Ohms 

RP 

parameter 

inductor effective resistance 

Ohms 


(primary) 


NP 

parameter 

no. of turn of primary 


NS 

parameter 

no. of turn of secondary 



OUTPUT 



V2 

variable 

output voltage 

Volts 


IP 

variable 

output current to source 

Amps 


PHI 

state variable 

flux 



VC 

state variable 

capacitor voltage 

Volts 


IS 

variable 

secondary current 

Amps 


VL 

variable 

inductor voltage 

Volts 



12 







EQUATIONS 



For IQ = 0, PHI > 0 : 

‘ IL I r-(RS+RC)/LS -1/NS] f IL ] [1/L RC/NS] f VI 

* = + 

. VC J L NS/C/LS 0 J L VC J [0 -1/C J 1 12 


‘ V2 ] NS*RC/LS 

IP J 0 



For IQ = 0, PHI = 0 ; 



13 



** LARGE SIGNAL MODEL OF FLYBACK (W/O LOAD) 
* POWER STAGE ONLY 

*#*****#*•*•■#■#■#••«•*•»•**•* *********»****' *•#••«■ •«••«• •«• ■«• ■«•*«■ 

** REVISED 12/7/85 

MACRO FILE NAME-MACROS 

DEFINE MACRO-FB 

MACRO INPUTS— 



C 

IP LS' NS NP 


RP 

RS RC 


OO 

VI 12 

* 

VI i 

INPUT VOLTAGE 

* 

12 ; 

INPUT CURRENT (FROM LOAD) 

# 

LP j 

PRIMARY INDUCTANCE 

*• 

LS i 

SECONDARY INDUCTANCE 

2 * 

NP j 

# OF TURN (PRIMARY) 


NS i 

# OF TURN (SECONDARY) 


QQ > 

SWITCHING FUNCTION 



QQ=1 (SW: ON) 



QQ=0 (SW: OFF ) 

MACRO OUTPUTS* 


PHI 11 IP IS 


VL 

VC V2 

* 

PHI i 

FLUX (STATE) 

CPACITOR VOLTAGE (STATE) 

* 

VC ; 


IP j 

PRIMARY CURRENT 


IS i 

SECONDARY CURRENT 


MACRO CODE 
MACRO STOP SORT 

I F < DABS ( GO FB — ) . LT. 1. E -5) THEN 
IF (PHIFB — . LF. 0. > THEN 
MACRO DER I VAT I VE, PH I FB— =0. 

GOTO <*+i-77 
ELSE 

TOFF ***** 1 

,A1 1=-(RS FB RC FB — )/LS FB — 

A12=-1/NS FB — 

A21=NS FB — /<C FB — *LS FB — ) 
A22=0. 

B11=0 

B 1 2=RC FB — /NS FB — 

B21=0. 

B22=-l/C FB — 

C11=NS FB — *RC FB— /LS FB — 

C 12=1 
C21=0 
C22=0 
D11=0 

D12=-RC FB — 

D21=0 

022*0 

* 

SE=1 

END IF 
ELSE 

*** TON *** 

All =— RP FB— /L.P FB — 

A 12=0 
A21=0 
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A 2 2=0 

bu-i/np r-n - - 

B 12=0 
B21=0 

B22=-l/C FB — 

C11=0 
Cl 2=1 

C21=NP FB—/LP FB — 
C22=0 
D 1 1=0 

D12=-RC FB — 

D21=0 
D22=0 
SE=0 
END IF 
** OUTPUT EQ. 

# 


* — - 


*\ 

V2 

• 

Cll 

ci 2 : 

! PHI! 

! Dll 

D12 I 

! VI 

I 

t 


* ! 


2= 


t 

» 

! ! + 

« 

1 

1 

1 
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C21 

C22 ! 

! VC ! 

! D21 

D22 ! 

! 12 



«. i i • « » » iii 

w i • i i i i iii 

V2 FB — =C11 *PHIFB — +C12*VC FB — 

S< + Dll #V1 FB — +D12#I2 FB — 

II FB — =C21*PHIFB — 

VL FB — =A11*PHIFB — +A12#VC FB <-Bll#Vl F6 — + 

& B12*I2 FB — 

IP FB — =1 1 FB — 

IS FB — =SE*NS FB — *PHIFB — /LS FB — 

**• STATE EQ. 

«• 



dPHI/dt 

! ! All 

A12 i 

! PHI ! 

! B 1 1 

B 12 ! i 



r: { 

I 

1 

I » 

) » 

<- i 

1 > 

1 \ 


' dVC/dt 1 

i A21 

A23 1 

: vc 

IB21 

B22 I 1 


3 

< 

1 

f 

1 

1 

i i 

» 1 

1 

3 

1 l 

» t 


«■ i 
* 


MACRO DERIVATIVE i PHIFB— — =VL FB — 

+++77 CONTINUE 

MACRO DERIVATIVE, VC FB — =A21*PHIFB — +A22*VC FB — + 
& B21-K-V1 FB — +B22*I2 FB — 

«• 

MACRO RESUME SORT 
END OF MACRO 
MODEL DESCRIPTION 
LOCATION= 22, FB 
END OF MODEL 
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<^-J 



VI 

variable 

input voltage 

Volts 

12 

variable 

input current from load 

Amps 

IQ 

variable 

switching function 




(if switch is on, IQ=1 ) 




(if switch is off, IQ=0) 


c 

parameter 

capacitance 

Farads 

L 

parameter 

inductance 

Henries 

RC 

parameter 

capacitor effective resistance 

Ohms 

RL 

parameter 

inductor effective resistance 

Ohms 

N 

parameter 

turns ratio(N=NS/NP) 



OUTPUT 


V2 

variable 

output voltage 

Volts 

IP 

variable 

output current to source 

Amps 

IL 

state variable 

inductor current 

Amps 

VC 

state variable 

capacitor voltage 

Volts 

IS 

variable 

switch current 

Amps 

VL 

variable 

inductor voltage 

Volts 
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EQUATIONS 


- (RL+RC)/L 


RC/L 

-i/c 


tc | ‘ 

o J 


For IQ “ 0, IL > 0 : 


- (RL+RC)/L 


For IQ = 0, IL « 0 


r 0 

L 0 


r 0 

L o 


17 



ft##*#**#***####*#**#*****##**####*#**#*##*# 

** LARGE SIGNAL MODEL OF FORWARD CONVERTER 
* POWER STAGE ONLY 

ft****************************************** 

** REVISED 12/20/85 
MACRO FILE NAME-MACROS 
DEFINE MACRO* FW 
MACRO INPUTS* 



C 

L RL RC 


IQ 

VI 12 


NS 

NP 

* 

IQ i 

SWITCHING FUNCTION 

# 


IQ=1 (SWITCH; ON) 

# 


IG=0 (SWITCH; OFF) 

** 

VI ; 

INPUT VOLTAGE 

* 

12 ; 

INPUT CURRENT (FROM LOAD) 

MACRO OUTPUTS-- 


IS 

IL 


ILL VL 


VC 

V2 IP 

* 

IL ; 

INDUCTOR CURRENT (STATE) 

* 

VC ; 

CAPACIIOR VOLTAGE (STATE) 

# 

IS ; 

SWITCH CURRENT 

# 

ILL, 

INDUCTOR CURRENT (DUMMY) 


VL ; 

TRANSFORMER VOLTAGE 

* 

V2 ; 

OUTPUT VOLTAGE 

* 

IP ; 

PRIMARY CURRENT (TO SOURCE) 


MACRO CODE 
MACRO STOP SORT 

XN=NS FW— /NP FW— 

IF <DABS( IQ FW — ). LT. 1. E-10HHEN 
IF ( IL FW— . LE. 0. ) THEN 
ILLFW — =0. 

MACRO DERIVATIVE, IL FW— =0. 

GOTO ♦•++•/ 7 
ELSE 

*#** TOFF ***** 

A 1 1=-(RL FW--+RC FW — >/l. FW— 

A 12= -1/1. FW- 
A21 = l/C FW— 

A22=0. 

B 1 1 =0 

B12=RC FW — /L FW— 

B21=0. 

B22=-l/C FW — 

Cl 1=RC FW— 

Cl 2=1 
C21=0 
C22=0 
D1 1=0 

D12=-RC FW— 

D21=0 

D22=0 

# 

IS FW — =0 
END IF 
ELSE 

*#* TON **# 

A1 1=— <RL FW — -t-RC FW — >/t FW — 
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A.t 2— ~1 /L FW - * 

A21= 1/C FW - - 
A22=0 

B 1 1 = XNl FW - - 
B12=RC FW — /L FW - - 
B21=0 

B22=— 1 /C FW — • 
CU=RC FW — 

C 12=1 
C21=1/XN 
C22=0 
D1 1=0 

D12=-RC FW — 

D21=0 

D22=0 

IS FW — = 1 L. FW — 

END IF 
** OUTPUT EQ. 

*• 


* 


*i V2 : 

*i ! = 

Cll 

C 12 i 

1 

1 

i IL 
» 

1 

b 

: Dii 

t 

i 

D12 ! 

t 

* 

: vi 

i 

i 

*: ip : 
* __ | 

C21 

C22 ! 
* 
I 

! VC 
1 
1 


! D21 
1 
1 

D22 : 
1 
1 

: 12 

i 

t 


* 


V2 FW — -Cl 1 *TL FW — +C12*VC FW — 

& + D11*V1 FW — +D12*I2 FW — 

IP FW — =C21*IL FW — 

VL FW — =<A11*IL FW — +A12*VC FW — +B11*V1 FW — + 
S< B 12*12 FW — >*L FW — 

** STATE EQ. 

* 


j dIL/dt 

j All 

A 12 i 
» 

! IL 

* 

{ B 1 1 

4- 1 

Bl? i 

I 

! VI ! 
) 1 

j dVC/dt 

! A 21 
1 
» 

A 22 ! 
1 
1 

1 VC 
1 
I 

! B 21 
1 
1 

B 22 ! 
1 
1 

! 12 ! 
1 I 

1 i 


MACRO DERIVATIVE, IL FW — =VL. FW— /L FW— 

-M-^-77 CONTINUE 

MACRO DERIVATIVE, VC FW — =A21*IL FW — +A22*VC FW — + 
& B21 *V 1 FW— +B22*I2 FW — 

& 

ILLFW — ~IL FW — 

IF< ILLFW — . I E. O. ) ILLFW— =0. 

MACRO RESUME SORT 
END OF MACRO 
MODEL DESCRIPTION 
LOCATIGN= 22, FW 
END OF MODEL 
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BATTERY DISCHARGER POWER STAGE 


BP 



INPUT 


VI 

variable 

input voltage 


Volts 

V2 

variable 

bus voltage 


Volts 

IQ 

variable 

switching function 
(if switch is on, 
(if switch is off. 

IQ-1) 

IQ-Q) 


L 

parameter 

inductance 

Henries 

RL 

parameter 

inductor effective resistance 

Ohms 


OUTPUT 


IL 

state variable 

inductor current 

Amps 


EQUATIONS 
For IQ = 1 : 

IL = ((1 + GAM) * VI - IL * RL - V2 ) / L 


For IQ = 0, IL > 0 : 

IL a (VI - IL * RL - V2) /L 

For IQ = 0, IL = 0 : 

IL = 0 
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■H it -ft- # 

* LARGE SIGNAL MODEL OF BATTERY DISCHARGER 

* POWER STAGE ONI Y 

* REVISED 3/4/86 
MACRO FILE NAME^MACROS 
DEFINE MACRO- BP 
MACRO INPUTS- 

L Rt. 

VI V2 V2L 

IG GAM EPS 

* IQ ; SWITCHING FUNCTION 

* IQ“1 (SWITCH; ON) 

* IQ=0 (SWITCH* OFF) 

* VI ; INPUT VTG FROM BATTERY 

* V2 i BUS VTG 
*• 

MACRO OUTPUTS- IL ILL Q 

* IL i INDUCTOR CURRENT (STATE VARIABLE) 

MACRO CODE 

MACRO STOP SORT 

q pp — -jq BP — 

IF(V2 BP — . GT. V2LBP — >G BP — =0 
IF ( Q BP — . LT. EPSBP — )THEN 
IF< IL BP — . LE. 0. )THEN 
ILLBP — -0. 

MACRO DERIVATIVE, IL BP~-*0. 

GOTO *■■*■*■ 77 

ELSE 

CC=1. 

END IF 
FL 8F 

CC=1. +GAMBP — 

END TF 

MACRO DERIVA1 1VE, IL BP — =(CC*V1 BP IL BP — *RL BP 

S< V2 BP — >/L BP — 

+++77 CONTINUE 

ILLBP — =-Jl BP — 

IF( ILLBP — . LE. O. ) ILLBP — =0. 

MACRO RESUME SORT 
END OF MACRO 
MODEL DESCRIPTION 
LOCATIQN= 22, BP 
END OF MODEL 
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LEAD/LAG COMPENSATOR 


DG 


V 


(K*VO-ER) 


GA (1+S/WZ) 
(1+S/WP) 


E 


INPUT 


VO 

variable 

input bus voltage 

Volts 

K 

parameter 

voltage dividing factor 


ER 

parameter 

op . amp . reference voltage 

Volts 

WZ 

parameter 

frequency of zero 

rad/sec 

WP 

parameter 

frequency of pole 

rad/sec 


OUTPUT 


VE 

variable 

ac component of control vtg. 

Volts 

XI 

state variable 

intermediate state variable 



EQUATIONS 


VE = [ XI + (K * VO - ER) * GA / WZ ] * WP 
X # 1 = GA * ( K * VO - ER ) 
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#**#***#***« 

*** COMPENSATOR (LEAD, LAG) 
tt*******#*^****#**^****#**#* 

MACRO FILE NAME-MACROS 
DEFINE MACRO=DG 
MACRO INPUTS- 

GA WZ WP 

K ER VO 

* 

* VO ; POWER STAGE OUTPUT VOLTAGE 

* WZ ; ZERO FREQUENCY < 2*P I #F ) 

* WP » POLE FREQUENCY ( 2*P I *F> 

* K » OUTPUT VOLTAGE DIVIDING RATIO 

* (IF VTG IVIDER NOT USED, K=l) 

* ER > REFERENCE VOLTAGE 

* 

MACRO OUTPUTS-X 1 VE 

* XI i DUMMY STATE 

* VE i OUTPUT VOLTAGE (TO PWM) 

MACRO CODE 

MACRO STOP SORT 

* 

TC1 = 1/W7. DG — 

TC2=1/WP DG — 

ERR=K DG—*VO DG — - ER DG — 

VE DG --(XI DG — +ERR*TC1*GA DG — )/TC2 
MACRO DERIVATIVE, XI DG--=GA DG— *ERR - VE DG 
# 

MACRO RESUME SORT 
END OF MACRO 
MODEL DESCRIPTION 
L0CATI0N=22, DG 
END OF MODEL 
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TWO-POLE ONE- ZERO COMPENSATOR 


MP 


V 

WM (1+S/WZ) 

V 

- O N 

(K*VO-ER) — 

E 


S (1+S/WP) 



INPUT 


VO 

variable 

input bus voltage 

Volts 

K 

parameter 

voltage dividing factor 

i 


ER 

parameter 

op. amp. reference voltage 

Volts 

WZ 

parameter 

frequency of zero 

rad/sec 

WP 

parameter 

frequency of pole 

rad/sec 


OUTPUT 


VE 

state variable 

ac component of control vtg. 

Volts 

XI 

state variable 

intermediate state variable 



EQUATIONS 


ERR= K * VO - ER 

XI = ERR - ( XI + ERR / WZ ) * WP 
VE = WM * ( XI + ERR / WZ ) * WP 
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*** COMPENSATOR <1 ZERO, 2 POLE) 

*4H**********##tt**#*******#***-* 

MACRO FILE NAME=MACROS 
DEFINE MACRO=MP 
MACRO INPUTS= 

WM WZ WP 

K ER VO 

* 

* VO { POWER STAGE OUTPUT VOLTAGE 

* WZ i ZERO FREQUENCY ( 2*P I *F> 

* WP {POLE FREQUENCY (2*PI*F> 

* K ; OUTPUT VOLTAGE DIVIDING RATIO 

* (IF VTG IVIDER NOT USED, K* 5 ! > 

* ER {REFERENCE VOLTAGE 

* 

MACRO OUTPUTS- XI VE 

* XI {DUMMY STATE 

* VE { OUTPUT VOLTAGE ( TO PWM) 

MACRO CODE 

MACRO STOP SORT 

* 

TC1 = 1/WZ MP-- 
TC2=1/WP MP — 

ERR=K MP- — *VO MP — - ER MP — 

MACRO DERIVATIVE, XI MP — =ERR - (XI MP — +ERR*TC1 > /TC2 
MACRO DERIVATIVE, VE MP — =QA- MP — *(X1 MP— +ERR*TC1) /TC2 

MACRO RESUME SORT ^ 

END OF MACRO 
MODEL DESCRIPTION 
LOCATION* 5 MP 
END OF MODEL 
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TWO-POLE TWO-ZERO COMPENSATOR 



(K*VO-ER). 


WM( 1+S/WZl ) ( 1+S/WZ2 ) 
S (1+S/WP) 


INPUT 


VO 

variable 

input bus voltage 

Volts 

K 

parameter 

voltage dividing factor 


ER 

parameter 

op. amp. reference voltage 

Volts 

WZ1, WZ2 

parameters 

frequency of zero 

rad/sec 

WP 

parameter 

frequency of pole 

rad/sec 


OUTPUT 


VE 

variable 

ac component of control vtg. 

Volts 

XI, X2 

state variables 

intermediate state variable 



EQUATIONS 

Z2 = WM * WP / WZ1 /WZ2 
ERR= K * VO - ER 
VE = X2 + Z2 * ERR 
XI = ERR * WM * WP 

X*2 = XI +Z2 * (WZ1 + WZ2 ) * ERR - WP * VE 
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*** COP ZENS ATOR (2 ZERO, 2 POLE) 
**#*************■**********■*■«•*•* 

MACRO FILE NAME -MACROS 
DEFINE MACROS Z 
MACRO INPUTS= 

WM WZ1 WZ2 WP 
K ER VO 

* 

* VO i POWER STAGF OUTPUT VOLTAGE 

* WZ1, WZ2 i ZERO FREGUENC I ES ( 2*P I *F ) 

* WP j POLE FREQUENCY (2*P I *F> 

* K i OUTPUT VOLTAGE DIVIDING RATIO 

* (IF VTG IVIDER NOT USED, K*1 > 

* ER i REFERENCE VOLTAGE 

* 

MACRO OUTPUTS- XI X2 VE 

* X 1 , X2 i DUMMY STATE 

* VE i OUTPUT VOLTAGE (TO PWM) 

MACRO CODE 

MACRO STOP SORT 
•* 

Z2-WM PZ— *WP PZ — / WZ 1 P Z — /WZ2P Z — * 

ERR=K P7 — *VO PZ ER PZ— 

V£ pz — -X2 PZ — +Z2*ERR 

MACRO DERIVATIVE, XI PZ --ERR *WM PZ — *WP PZ — 

MACRO DERIVATIVE, X2 PZ — =X1 PZ — +Z2*(WZ1PZ — +WZ2PZ — )*ERR 
& -WP PZ— *VE PZ — 

#• 

MACRO RESUME SORT 
END OF MACRO 
MODEL DESCRIPTION 
LOCATION-22, PZ 
END OF MODEL 
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S.C.M CURRENT LOOP 


SM 



INPUT 


mm 

variable 

input bus voltage 

Volts 


parameter 

transformer turns ratio 


c 

parameter 

current- loop capacitance 

Farads 

R 

parameter 

current-loop resistance 

Ohms 


OUTPUT 


VI 

state variable 

current- loop output voltage 

Volts 


EQUATIONS 


VI = NV * VL / ( C * R ) 


28 












********************************#.*.#■# 

*** CURRENT LOOP (S C M> MODULE 
*#******«***************«*#*-*#«•*•»#•*#*### 

* 

MACRO FILE NAME-MACROS 
DEFINE MACRO^SM 
MACRO INPUTS= 

VL NV C R 

# 

* VL i TRNASFORMER VOLTAGE (FROM POWER STAGE) 

* NV i TRANSFORMER TURNS RATIO 

* 

MACRO OUTPUTS=VI 
* 

* VI i CURRENT LOOP OUTPUT VOLTAGE 

* (TO PWM BLOCK) 

MACRO CODE 
MACRO STOP SORT 

ESTATE EQUATION 

MACRO DERIVATIVE, VI SM — =NV SM — *VL SM — / ‘ 

& (C SM — *R SM— ) 

*• 

MACRO RESUME SORT 
END OF MACRO 
MODEL DESCRIPTION 
LOCATION* 1 22, SM 
END OF MODEL 
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CURRENT. INJECTION CONTROL 



INPUT 


IL 

variable 

inductor current 
(if flyback, flux) 

! 

Amps 

NI 

parameter 

current transformer tur;ns ratio 


NP 

parameter 

if flyback, primary # of turns 
otherwise NP=1 


LP 

parameter 

if flyback, primary inductance 
otherwise LP=1 

Henries 

RW 

parameter 

current loop resistance 

Ohms 


OUTPUT 


VI 

variable 

current loop output voltage 

Volts 


EQUATIONS 

For FLYBACK CONVERTER : 

VI = NP / LP * IL * RW / NI 

For OTHER CONVERTERS : 

VI = IL * RW / NI 
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*■** CURRENT LOOP <C I C) MODULE 
*****#**********##*##*###*#***#.*###*. 

MACRO FILE NAME^MACROS 
DEFINE MACRO=CC 
MACRO INPUTS= 

IL NI RW 

NP LP 

* 

* IL i INDUCTOR CURRENT (FROM POWER STAGE) 

* (FLUX IF BUCK/BOOST) 

* NI ; CURRENT TRANSFORMER TURNS RATIO 

* NP i IF FLYBACK : POWER STAGE PRIMARY # OF TURNS 

* OTHERWISE NP=1 

* LP i IF FLYBACK ; POWER TRANSFORMER PRIMARY INDUCTANCE 

* OTHERWISE LP=i 

* 

MACRO OUTPUTS=VI 
* 

* VI i CIC OUTPUT CONTROL VOLTAGE 

* 

MACRO CODE 
MACRO STOP SORT 
■* 

VI CC — ~NP CC- /LP CC--* 

& IL CC — *RW CC- -*/NI CC — 

* 

MACRO RESUME SORT 
END OF MACRO 
MODEL DESCRIPTION 
L0CATI0N=22j CC 
END OF MODEL 
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PWM( CONSTANT FREQUENCY CONTROL) 


WM 



IQ 
- — ► 


INPUT 


VE 

variable 

voltage -loop error voltage 

Volts 

VI 

variable 

current-loop error voltage 

Volts 

TI 

parameter 

switching period 

Seconds 

VP 

parameter 

amplitude of external ramp 

Volts 

VQ 

parameter 

thershold voltage 

Volts 

ER 

parameter 

reference voltage of op. amp. 

Volts 


OUTPUT 


VR 

variable 

reference voltage 

Volts 

VC 

variable 

total control voltage 

Volts 

IQ 

variable 

switching function 



EQUATIONS 

VR = VP * (TN - N) + VQ 
VC = - VE - VI 

IF VR > VC , IQ = 0 
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********«###*##****#»«# »*#•»«•*« tHf-tt# 

** P. W. M (CONSTANT FREQUENCY ) *■*£• 

•**•#•# -fi- #***•###***•#* tf###**#**'**#*** 


MACRO FILE NAME- MAC R OS 
DEFINE MACRO" WM 


MACRO INPUTS- 


TI 

VP 

VQ 

CIC 

SCM 


VE 

VI 

EPS 

vex 

VCN 

I LX 

DMN 

DMX 



* 


* VE ; AMPLIFIED ERROR VTG ( FROM COMPENSATOR) 

* VI ; INPUT VOLTAGE FROM CURRENT FEEDBACK MODULE 

* CIC* IF CIC, CIC=1 

* OTHERWISE CIC=-0 

* SCMi IF SCM* SCM=1 

* OTHERWISE SCM=0 

* TI * SWICHING INTERVAL 

* VP i AMPLITUDE OF EXTERNAL RAMP 

* ( IF EXT. RAMP NOT USED* VP-O) 

* VQ ; THRESHOLD VOLTAGE 

* ER i REFERNCE VOLTAGE OF OP. AMP 

* 


MACRO OUTPUTS^ 



VR 

VC IQ TN 

& 

VS 


# 

VR i 

EXTERNAL RAMP VOLTAGE 


VC ; 

TOIAL CONTROL VOLTAGE 

# 

IQ i 

SWITCHING FUNCTION 

* 


(IP SWITCH=ONi IQ-1) 

* 


(IF SWITCH=OFF* IQ^-O) 


MACRO CODE 
MACRO STOP SORT 

*** RAMP GENERATION 
*• 

TN WM — ~ < T I ME +T I WM — >/TI WM — 

NP=N 

N=IDINT < TN WM — ) 

VR WM — ~VP WM— *(TN WM- N)+VO WM — 

VC WM — =. -VE WM CICWM — *VS WM 

& SOMWM TfVI WM — 

* OP-AMP SATURATION 

IF< VC WM — . GT. VCXWM — >VC WM — =VCXWM— 

IF ( VC WM — . LT. VCNWM — >VC WM — “VCNWM — 

* COMPARATOR 

IF < IQ WM— . GT. EPS) THEN 
VS WM — ~ V3 WM — 

IF < VR WM— . GT. VC WM — ) THEN 
IQ WM— -O 
END IF 
ELSE 

VS WM — -0. 

END IF 

IF ( NP. NE. N) IQ WM — =1 
^PROTECTION 

D=TN WM' N 

IF ( D. LT . DMMWM — ) 1 0 WM — = i 
IF ( D. GT. DMXWM— ) TQ WM — ~0 
IF ( IL WM— . Cl. ILXWM — )IQ WM — =0 

MACRO RESUME SORT 
END OF MACRO 
MODEL DESCRIPTION 
LOCATION= 22 * WM 
END OF MODEL 
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PWM( CONSTANT OFF-TIME CONTROL) 




INPUT 


VE 

variable 

voltage-loop error voltage 

Volts 

VI 

variable 

current-loop error voltage 

Volts 

TOF 

parameter 

constant off-time interval 

Seconds 

SLP 

parameter 

slope of external ramp 

Volts/sec 

VTH 

parameter 

threshold voltage 

Volts 

ER 

parameter 

reference voltage of op. amp. 

Volts 


OUTPUT 


Q 

variable 

switching function 


VR 

variable 

reference voltage 

Volts 

VC 

variable 

total control voltage 

Volts 


EQUATIONS 
VC = - VE - VI 

VR = V(RAMP) + VTH 
IF VR > VC, Q = 0 
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# P W M (CONSTANT OFF TIME CQNTROl > 

MACRO FILE NAME-MACROS 
DEFINE MACRO- FF 
MACRO INPUTS- 

VE VI IL 

TOF EPS VTH SLP ILM 

ER CIC SCM 

# 

* VE ; AMPLIFIED ERROR VTG(FROM COMPENSATOR) 

* Vis INPUT VTG FROM CURRENT FEEDBACK MODULE 

* CICi IF CIC* CIC=1 

* OTHERWISE CIC=0 

*• SCMi IF SCM, SCM=1 

* OTHERWISE SCM-O 

* TOF i CONSTANT OFF TIME INTERVAL 

* SLP ; SLOP OF EXTERNAL RAMP 

* ( IF EXT. RAMP NOT USED, SLP=0) 

* VTH; THRESHOLD VOLTAGE 

* ER ; REFERNCE VOLTAGE OF OP. AMP 

* 

MACRO OUTPUTS- 

VR VC Q 

* VR ; EXTERNAL RAMP VOLTAGE PLUS THRESHOLD VTG 

* VC ; TOTAL CONTROL VOLTAGE 

* G ; SWITCHING FUNCTION 

* (IF SWITCH^ON, G~1 ) 

* (IF SWITCH=QFF, Q—0) 

* TC ; OFF TIME INSTANT 

* TR ; ON TIME INSTANT 

MACRO CODE 

MACRO STOP SORT 

fi- 

VC FF— -- -VE FF~- — CICFF — #VI FF 

& SCMFF — *VI FF — 

IF (TIME. GT. EPS) GOTO +++1 1 
T-TIME 
Q FE-~- 1 
GOTO +++22 
+++11 CONTINUE 
*• 

IF(Q FF — . LT. EPSFF — > THEN 

«■ 

.IF ( (TIME-TC). GE. TOFFF — ) THEN 
TR--TIME 
G FF— =1 
END IF 

ELSE 

VR FF — =SLPFF— * (TIME -TR > +VTHFF — 

IF ( VR FF — . GE. VC FF — >THEN 
Q FF---0 
TC-TIME 
END IF 

END IF 

+++22 CONTINUE 

I F ( I L FF - -. GE. ILMFF — ) G FF — -O 
MACRO RESUME SORT 
END OF MACRO 
MODEL DESCRIPTION 
LOCATION^ 22, FF 
END OF MODEL 
# 
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PWM( CONSTANT ON-TIME CONTROL) 



NN 


INPUT 


VE 

variable 

voltage loop error voltage 

Volts 

VI 

variable 

current loop error voltage 

Volts 

TOF 

parameter 

constant off-time interval 

Seconds 

SLP 

parameter 

slope of external ramp 

Volts/sec 

VTH 

parameter 

threshold voltage 

Volts 

ER 

parameter 

refernce voltage of op. amp. 

Volts 


OUTPUT 


Q 

variable 

switching function 


VR 

variable 

reference voltage 

Volts 

VC 

variable 

total control voltage 

Volts 


EQUATIONS 
VC = - VE - VI 

VR = V(RAMP) + VTH 
IF VR < VC, Q = 1 
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*##***#*##**#«** tt-*****-* 

* P W M (CONSTANT ON TIME OONTROl ) 

MACRO FILE NAME- MACROS 
DEFINE MACRO" NN 
MACRO INPUTS- 

VE VI IL 

TON EPS VTH SLP 1LM 

ER CIC SCM 

*• 

* VE i AMPLIFIED ERROR VTG(FROM COMPENSATOR) 

* VI i INPUT VTG FROM CURRENT FEEDBACK MODULE 

* CIC; IF CIC* CIC=1 

* OTHERWISE CIC-0 

* SCM; IF SCM, SCM=1 

* OTHERWISE SCM=0 

* TON ; CONSTANT ON TIME INTERVAL 

* SLP ; SLOP OF EXTERNAL RAMP 

* ( IF EXT. RAMP NOT USED. SLP=0) 

* VTH; THRESHOLD VOLTAGE 

* ER ; REFERNCE VOLTAGE OF OP. AMP 

* 

MACRO OUTPUTS^ 

VR VC Q 

* 

* VR ; EXTERNAL RAMP VOLTAGE PLUS THRESHOLD VTG 

* VC ; TOTAL CONTROL VOLTAGE 

* Q ; SWITCHING FUNCTION 

* (IF 5WITCH=0N* G=!> 

* (IF SWITCH=QNN» Q-O) 

* TC ; ON TIME INSTANT 

* TR ; OFF TIME INSTANT 

MACRO CODE 

MACRO STOP SORT 

*• 

VC NN — = -VE NN— Cl CNN — *VI NN 

& SCMNN — *VI NN — 

IF(TIME. GT. EPS)GOTO +++11 
T=TIME 
Q NN- "" I 
GOTO +++22 
+++11 CONTINUE 

ft 

IF(Q NN — . GT. EPSNN — )THEN 

* 

IF( (TIME-TC). GE. TONNN — ) THEN 
TR-TIME 
Q NN— =0 
END IF 

ELSE 

VR NN — =SLPNN — * ( T I ME -TR ) +VTHNN — 

IF ( VR NN — . LE. VC NN — >THEN 
Q NN— =1 
TC-TIME 
END IF 
END IF 

+++22 CONTINUE 

IF( IL NN — . GE. ILMNN — )Q NN— =0 

* 

MACRO RESUME SORT 
END OF MACRO 
MODEL DESCRIPTION 
LOCATION= 22, NN 
END OF MODEL 
x 
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INPUT 


VI 

variable 

input bus voltage 

Volts 

RA 

parameter 

resistance 

Ohms 

RB 

parameter 

resistance 

Ohms 

L 

parameter 

inductance 

Henries 

C 

parameter 

capacitance 

Farads 

TC 

parameter 

time for step change 

seconds 


OUTPUT 


IL 

state variable 

inductor current 

Amps 

VC 

state variable 

capacitor voltage 

Volts 

R 

variable 

resistance 

Ohms 


EQUATIONS 

For TIME < TC : R = RA 

For TIME > TC : R = RB 

IL = ( VI - VC ) / L 

VC = ( IL - VC / R ) / C 
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it-***#********#**##**#***#***##*# 

******** LOAD MODEL ********** 
***************************** 

MACRO FILE NAME- MACROS 
DEFINE MACR0-L.0 
MACRO INPUTS^ VI RA RB 

L C TC 

* VI ; INPUT VOLTAGE (FROM POWER STAGE) 

* RA i RESISTOR VALUE BEFORE STEP CHANGE 

* RB i RESISTOR VALUE AFTER STEP CHANGE 

* TC J TIME FOR STEP CHANGE 

MACRO OUTPUTS- IL VC R 

MACRO CODE 
MACRO STOP SORT 

R LO — -RA LO~ 

IF (TIME. GT. TC LO — >R LO— -RB LO— 

*ST ATE EQUATION 

MACRO DERIVATIVE, I L LO — -(VI LO — - VC LO — ) /L LO — 

MACRO DERIVATIVE, VC LO— =<IL LO VC LO — /R LO — ) /C LO 

* 

MACRO RESUME SORT 
END OF MACRO 
MODEL DESCRIPTION 
L0CATI0N=25, LO 
END OF MODEL 
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CONSTANT POWER LOAD 




INPUT 


VL 

variable 

input bus voltage 

Volts 

PW1 , PW2 
PWO 

parameters 

constant power values 

( 

Watts 

PC 

parameter 

if PC-1, const, power 
if PC=0, time varying power 
if PC=2, time varying power 


VR 

parameter 

minimum voltage to 
maintain const. power 

Volts 

SW 

parameter 

slope of time varying power 

Watts/sec 


OUTPUT 


IL 

variable 

load 

current 

Amps 

PW 

variable 

load 

power 

Watts 


EQUATIONS 


For PC = 1 : 

PW = PW1 ( IF TIME < TC ) 
' PW = PW2 ( IF TIME > TC ) 


For PC = 2 : 

PW = PW1 - ( SW * TIME ) 
PW = PW2 ( IF PW < PW2 ) 


For PC = 0 : 

PW = SW * TIME + PWO 
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MACRO FILE NAME *= MACROS 

ft •»■&•}!• -ft* 

DEFINE MACRO =- El 
ft LOAD MODEL 

* INPUT IS BUS VOLTAGE <VL) AND POWER <PW> OUTPUT IS I (LOAD) 

ft 

ftftftft FOR STEP POWER CHANGE s PC PT = 1. 

**** FOR RAM P CHANGE OF POWER i PC PT -0. OR 2. 

ft PC PT ~ 0 : PW = SW ft TIME + PWO 

* PC PT - 2 : PW - PWi - (SW ft TIME) 


MACRO 

INPUTS = 

Ut 

PO 

PWI 

PW? 

PC 



VR 

TC 

SW 



MACRO 

OUTPUTS »■ 

IL 

PW 

SL 




MACRO CODE 
MACRO STOP SORT 


IF ( PC PT — . EQ. 1. ) 
IF ( PC PT ~ . EQ. 2. ) 
P W P T — SW PT — ftTI ME 

IL PT — - PW PT — / VL 
GOTO + + +33 


GO TO +++23 
GO TO +++43 
+ PWOPT — 

PT — 


•ft STEP CHANGE OF POWER AT TIME = TC PT 
+++23 CONTINUE 

PW PT — -PW1PT — 

IF (TIME. GE. TC PT — )PW PT — -PW2PT — 

IL PT --- =• PW PT — / VL PT — 

+++33 SL PT — -• PW PT — /VR PT — /VR P7 — 

IF (VL PT — . LT. VR PT— > IL PT— ~ SL PT — *VL PT — 
IF ( IL PT— - . LT. O. > IL PT — » O. 

GO TO I- + +53 


* RAMP CHANGE OF POWER FROM PW1PT TO PW2PT WITH THE SLOPE 
ft OF SW PT C.W3 PER SECOND 

+++4.3 PW PT - - « PW j FT ( SW PT — ft T IME ) 

IF ( PW PT— . LT. PW2PT — ) PW PT — — P WOP T — 

IL PT — « PW PT — / VL. PT — 

SL PT— - PW PT- — / VR PT— / VR FT — 

IF ( IL PT— . LT. O. > JL PI — o. 

ft DUMMY STATEMENT 
+ + +53 PW PT— -- PW PI — 

MACRO RESUME SORT 
END OF MACRO 

ftftftftftftftftftftftftftftftftft#ftftftftftftftftftftftftftftftftftftftftftftft'ftftftftftftft,ii.fi-.*i-.ft-»ft-K- 

MODEL BE SCR I P T ION 
L.OC AT 1 0N=4 2 . PT 
END OF MODEL 
Print 
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SOLAR ARRAY [AR] 


AR 



INPUT 


Physical 

Quantity 

Name. ....... 


Description 

units 

IX 

variable 

Bus current 

amps 

LL 

parameter 

Illumination level ’ 


R 

parameter 

Filter resistance 

ohms 

L 

parameter 

Filter inductance 

henries 

C 

parameter 

Filter capacitance 

farads 

TA 

parameter 

Ambient temperature 

°K 

NP 

parameter 

No. of parallel arrays 


NS 

parameter 

No. of serial cell 


FC, FV 

parameter 

Temperature coeffieients 



OUTPUT 


* 

VB 

variable 

Bus voltage 

volts 

VO 

variable 

Solar array output voltage 

volts 

11* 

variable 

Filter inductor current 

volts 


* state variable 
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MACRO FILE NAME = MACROS 
DEFINE MACRO « AR 

*ft*ftftft***ftftftftftftftft*ftft*ft***ftft*ftft 

ftft SOLAR ARRAY MODEL **-**#* 

«•!- = = ====: = > VP 

#=======?==:«:< JX 

•ft -ft* 

MACRO INPUTS ~ IX L.L LL..S 

R L. C 1 A 

FC FV NP NS 

ft IX « INPUT CURRENT 
ft L.L - ILLUMINA'I JON LEVEL 
ft LLS - SLOPE OF ILLUMINATION CHANGE 
ft R» L, C, : CABLE IMPEDANCE 

* TA = AMBIENT TEMPERATURE 

* FC, FV : TEMPERATURE COEFFICIENTS 
ft NP = NO. OF PARALLEL ARRAYS 

* NS = NO. OF SERIAL CELLS 

* R, L, C» : CABLE IMPEDANCE 

f 

MACRO OUTPUTS - VO IL. VB 

* VO - SOLAR ARRAY OUTPUT VOLTAGE 

* VB -• BUS VOLTAGE 

MACRO CODE 
MACRO STOP SORT 

ftftft SOLAR CELL PARAMETERS **** 

ft- 
's- R5 = INTERNAL. SERIES RESISTANCE 

* RSH“ INTERNAL SHUNT RESISTANCE 

* XIO= REVERSE SATURATION CURRENT 

* IN - NOMINAL TEMPERATURE 

* Q ~ ELECTRON CHARGE 

* XK = BOLTZMANN CONSTANT 

* XIG= LIGHT-GENERATED CURRENT 

* VOC^ OPEN CIRCUIT VOLTAGE 

ft 

RS = .42 
RSH =• 2 DO. 

X IG = . H115 
XIO = 4. IS69E-11 
A ~ . 967 
TN ~ 30 J. 

Q -- 1. 60PE-19 
XK * 1 . 38 1 E--23 
VOC “ . 5012 

* XKO = 0 / ( XK ft A ft TN > 

XK0=39. S 

•ft 

* ILLUMINATION CHANGE WITH A LINEAR SLOPE OF LLS 

FILL = LL AR — r LL.SAR— - * 1 IME 

ft EFFECT OF ILLUMINATION CHANGE 


[ AR . MOD ] continued 
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Ir 


X1G ~ X I G ft FILL 


Cl - < 1. + RS / RSH > 

C2 = NP AR — / < NS AR — ft RSH ) 

C3 = -NP AR — * X 16 
A1 = XKO / NS AR — 

A2 = XKO * RS / NP AR- 
IF ( TIME . NE. O. ) 6010 +++10 

* INITIAL GUESS FOR SOLAR ARRAY OUTPUT VOLTAGE 

VO AR — » VB AR— 

ft##*#*#**#******####*#***##**#*#****#***######*###** 

### NEWTON ITERATION **#**#**###****#**#*#*##**# 

ft***#######***#*####*#####*****##***#***#*########*# 

f++10 CONTINUE 

VOP * VO AR — 

FV = Cl ft IL AR— * C2 * VO AR — + C3 + NP AR — 

& * XIO ft DEXP < A1 * VO AR — + A2 ft IL AR — ) 

DFV - C2 + Ai ft NP AR — ft XIO ft DEXP < A1 ft VO AR — 

S< + AL' ft IL AR— ) 

VO AR— =• VO AR- FV / DFV 

11 = < VO AR — - VOP ) / VO AR- 
IF ( DADS ( 11 > , L.E. 1. E-4 ) GOTO ++ t-20 

GO TO +++.10 

###########*##*#*•#•##**#***#**#*•#**#**##*#*#***#**##*#*## 
ft TEMPERATURE CORRECTION FOR S. A. I-V CURVE ft 

ft DEL ( I ) = FC ft DEL < T ) ft 

ft DEL < V ) = ( PV + FC ft RS ) ft DEL < T > ft 

*######***#######**##*#*##*#ft#***#**###**####*#*####*### 

+*+20 DELI ~ FC AR— ft < TA AR- TN ) 

DELV = ( FV AR — + FC AT! — ft RS ) ft < TA AR— - TN ) 

ILT = IL. AR — + DELI 
VO AR— * VO AR— + DELV 

# 

VL = VO AR — - I LI * R AR VB AR — 

MACRO DERIVATIVES, IL. AR— = VL / L. AR— 

MACRO DERIVATIVES, VB AR — = < IL AR- IX, AR— ) / C AR- 

MACRO RESUME SORT 
END OF MACRO 

**ft#**###*#*###*#*####ftft*#*****#*##*************####**#*#* 

MODEL DESCRIPTION 
LOCATION = 22, AR 
END OF MODEL 
PRINT 


[ AR . MOD ] 
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Physical 

Quantity 

Name 


variable 

variable 

parameter 

parameter 

parameter 

parameter 

parameter 

parameter 

parameter 

parameter 

parameter 


Description 

units 

Load current 

amps 

Shunt current 

amps 

f 

Delay time between switchings 

seconds 

Filter resistance 

ohms 

Filter inductance 

henries 

Filter capacitance 

farads 

Filter resistance 

ohms 

Filter capacitance 

No. of parallel arrays 

No. of serial cells 

No. of parallel arrays switched at a 
time 

farads 


OUTPUT 


VB* 

variable 

Bus voltage 

volts 

VO 

variable 

Solar array output voltage 

volts 

VC1* 

variable 

Filter capacitor voltage 

volts 

IS 

variable 

Solar array output current 

amps 

IL* 

variable 

Filter inductor current 

amps 

NPP 

variable 

No. of parallel arrays on bus 


* state variable 
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OUTPUT 


Physical 

Quantity 

Name 


i 

Description 

units 

VB* 

variable 

Bus voltage 

■1 

VOl 

varialbe 

Output voltage to upper array 

ns 

V02 

variable 

Output voltage of lower array 

volts 

IL1* 

variable 

Upper array output current 

amps 

12 

variable 

Lower array output current 

amps 


* state variable 
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ORIGINAL PAGE- IS 
OF POOR QUALITY 


MACRO FILE NAME * MACROS 
DEFINE MACRO - AS 

«-^*****^*«**** &»«•«■ •»•■»■«• *##*##•£*•* -iHfr**##* tt-S-H-S-S- 

** SOLAR ARRAY SWITCHING UNIT MODEL «•*•* 

KSPBBWSS:^ IX 


#*•& 


MACRO INPUTS ~ jx IH DL 

R L C R i C i 

NP NS NNP 

* IX a LOAD CURRENT 

* IH - SHUNT CURRENT 

* DL » DELAY 1 1 ME BETWEEN SWITCHINGS 

* R> l.» C> : CABLE IMPEDANCE 

* Rl,Cl> : FILTER 

* NP » NO. OF PARALLEL ARRAYS 

* NS - NO. OF SERIES CELL NO. 

* NNP = NO. OF PARALLEL ARRAYS SWITCHED AT A TIME 


MACRO OUTPUTS -- VO IL VB NP P 

vei io 

* VO a sot AR ARRAY OUT VOLTAGE 

* VB ~ BUS VOLTAGE 

* NP.P = NO. OF PARALLEL ARRAYS 

* 10 a SOLAR ARRAY OUTPUT CURRENT 


MACRO CODE 
MACRO STOP SORT 


■*** SOLAR CELL PARAMETERS **** 

* RS « INTERNAL SERIFS RESISTANCE 

* RSH-“ INTERNAL SHUNT RESISTANCE 

* XI-O* REVERSE SATURATION CURRENT 

* TE * TEMPERATURE 

* Q ~ ELECTRON CHARGE 

* XR » BOLTZMANN CONSTANT 

v- XI G» SHORT CIRCUIT CURRENT 

* VGCa OPEN CIRCUIT VOLTAGE 

RS--. 42 
RSH=250. 

X I Q— . 14.1 15 
XI 0=4. 1 869E-.il 
A--. 969 
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TE~30i . 

Q“i. 602E - .1? 

XK-1 . 381E-23 
VOC=. bb!2 

* XKG-Q/<XKftAftXE) 

XKQ-3?. 8 

ftftftftft^frftftftftftftftftftftftftftftftftftftftftftftftftftftftftftftftftftftftftftft *• -ft •&«••»?•«• R-5H5- -ft *■»■# •«•■»•■«••!«• -li- 
ft SOLAR ARRAY SWITCH I WO LOGIC ft 

ftf.-K"S$"K ftK--ft-*-ft-fr*-ft-ft*-*-ftft*fcfti>ftftft-ft#-ftft*ftft**ft**ftftftftft-Rft#*ftftftftftftft*#ftftftftftftftft 

IF (TIME. GT. 0. 5 GOTO +Hi 0 
NPPAS — =•• NP AS — 

F ! I Mi.— 0. 

CONTINUE 

IF < IH AS — . ©E. 5. .AND. NPPAS— . GT. 64 ) GOTO + >+20 
IF ( IH AS— . LT. 2. .AND. NPPAS-- .IT. 324 > GOTO +4- >30 
GOTO +M40 

IF < TIME . LT. FTIMEMDL AS— > G07 0 +++40 

NPPAS — =•• NPPAS- NNPAS — 

FT I ME «- TIME 
GOTO + + >40 

IF < TIME . I T. FTIME+DL AS — > GOTO ftft+40 
NPPAS— f NPPAS — + NNPAS — * 7 

FT IMF " TIME 
GOTO ftft+40 
CONTINUE 


ft ILLUMINATION CHANGE 
LL AS-- = 1. 

ft 

XIG - XIG ft LL. AS — 

Cl = •< i. + R3/RGH ) 

C2= NPPAS — / < NS AS — ft RSH ) 

C3= -NPPAS— ft XIG 
At- XKO / NS AS — 

A 2= XKO ft RS / NPPAS— 

IF ( TIME .ME. 0. ) GOTO ft ft ►50 

ft 

ft INITIAL GUESS 

VO AS — -VB AS— 

k ftft NEWTON ITERATION *ftft*ftftftft*ft*ftft*ftft 

ftft-R ftftR^ftftftftftftft-R ft-ftft-ftft-ftft-ft'ft-ft-ftftft- 

Eft +50 CONTINUE 

VOP-VO AS — 

ft I - V EQUATION 

ft CAPACITOR ADDED 

IF ( TIME. NF. 0. ) GOTO +++60 
10 AS— * 30. 5 

+++60 FV = Cl ft ID AS — > C2 ft VO AS — ft C3 + NPPAS— ft XI 0 ft 

& DEXP < At * VO AS— + A 2 ft JO AS— > 

DFV - C? ft A.1 ft NPPAS— ft XI 0 ft 
& DEXP < At ft VO AS — ft A2 ft 10 AS — ) 

VO AS— - VO AS — - FV / DFV 

72 - < VO AS — - VOP > / VO AS — 


■"ft Ft 0 

+++20 

+++30 

>++40 
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ORSGIMAL PAGE 1 IS 
OF POOR QUALITY 


IF ( DABS ( ZZ ) . LE. l.E-4 > OOTQ + M- 70 
GO TO M i-50 

+•+•+70 VL » VO AS- - - IL AS — #R AS — - VB AS— 

MACRO DERIVATIVES, VC 1 AS— =' (10 AS IL AS — >/Cl AS — 

MACRO DERIVATIVES, IL AS— = VL / L AS — 

MACRO DERIVATIVES, VB AS— = <IL AS— - IX AS— > /C AS — 
10 AS — - (VO AS — - VC I AS— ) / Ri AS— + IL AS— 
IF (VO AS— . LE. 0. ) VO AS— & 0. 

IF (in AS— . LE. 0. ) 10 AS— ~ 0. 

MACRO RESUME SORT 
END OF MACRO 

«•!+•«• •«••£*■» •ft**#**#***#*# ■«• 

MODEL DESCR I PT I DM 
LOCATION - 22, AS 
END OF MODEL 
PRINT 
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Physical 
Quantity 

Name I I Description I Units 


VB 

variable 

input bus voltage 

Volts 

VR 

parameter 

reference voltage 

Volts 

K1 

parameter 

bus voltage divider 
ratio 

Volts 

Al, ZO, PO 

parameters 

As block parameters 
(See figure) 


LI 

parameter 

OP-Arap saturation voltage 

Volts 

Z1,P1 

parameters 

G1 block parameters 
(See figure) 


A2 , Z2 , P2 

parameters 

G2 block parameters 
(See figure) 


VF 

parameter 

junction voltage drop of 
darlington circuit 

Volts 

K2 

parameter 

gain of control voltage to 
base current 


K3 

'1 

parameter 

No. of shunt elements 


( 

L3 

parameter 

Ic vs V(be) characteristic 

Volts 

Cl, C2 

1 

parameters 

shunt transistor circuit 
coefficients 


IMA 

parameter 

maximum shunt current 

Amps 

OUTPUT 

Vol 

variable 

output voltage of error 
amplifier 

Volts 

XI , X2 , S2 

variables * 

dummy state variables 


IH 

variable 

shunt current 

Amps 


* state variable 
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bus 
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* SHUNT ELEMENTS 



ORIGINAL PAGE IS 
OF POOR QUALITY 


MACRO FILE NAME = MACROS 

*****##*###*##*#*#**#***##*##***###■& ##**#*#*-*##***■»• 

# SHUNT REGULATOR MODEL < TYPE 1) WITH * 

* TRANSFER FUNCTION * 


DEFINE MACRO 

« HA 




MACRO INPUTS 

- VB 

VR 

K1 


At 

70 

PO 

l 1 


Zi 

P 1 

A2 

22 

P.2 

VF 

IMA 

K2 

K3 

1.3 

Cl 


* VB * BUS VOLTAGE 

* VR - REFERENCE VOLTAGE 

* K1 = BUS VOLTAGE DIVIDER RATIO 

* IMA = SHUNT MAXIMUM CURRENT 

MACRO OUTPUTS ~ VOl XI 

IH SR X? 

* VOl = OUTPUT VOLTAGE OF ERROR AMPLIFIER 

* XI, S2, X2 = DUMMY STATE VARIABLES 
*• IH ~ SHUNT CURRENT 

MACRO CODE 
MACRO STOP SORT 

* , — . — 

s- 

* THIS SECTION SIMULATES THE SHUNT ERROR AMPLIFIER 

* REFERENCE VOLTAGE IS VR HA— 

K- 

* « AS BLOCK » 

* ERROR AMPLIFIER MODEL 

GAI = A1 HA — 

SI = K1 HA — * VB HA — - VR HA- 

SS = GAI * SI + XI HA — 

MACRO DERIVATIVE, XI HA. — = GAI * SI * ZO HA — - 32 * PO HA — 

DD — S2 + VR HA — 

k- « LI BLOCK » 

* OP AMP OUTPUT VOLTAGE IS LIMITED TO A VALUE BETWEEN Li AND 
■* (VB - LI) 

BB = VB HA — - LI HA — 

CC - 11 HA — 

IF ( DD . L.T. CC > DD = CC 
IF ( DD . GT. BB ) DD = BB 
VOl HA— - DP 

* <:< Qi liL DCK » 

* MODEL OF NETWORK AT THE OUTPUT OF ERROR AMPLIFIER 

MACRO DERIVATIVE, SS HA — = Z1 HA — * VOIHA — - PI HA — * SS HA — 
W1 = SS HA— 

* « G2 BLOCK » 

GAI = AS HA- 
SS = VB HA — 
g4 =: GAI # S3 + X2 HA — - 

MACRO DERIVATIVE, ' X2 HA — = GAI * S3 * Z 2 HA 34 * P2 HA — 
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VV2 = 34 

>'* 

VGC -• VVJ <- W2 

Hr 

* MODEL OF TRANSISTOR DRIVER CIRCUIT 

VG = VB HA — - (VOC + VF HA— > 

* « L2 BLOCK » 

IF (VG. LT. 0. ) VQ - 0. 

* MODEL OF SHUNT ELEMENTS 

* « K2, K3 BLOCKS » 

FIB * VQ * K2 HA — * K3 HA — 

* « L3 BLOCK » 

IF < FIB . LI . L3 HA — ) FIB = 0. 

<- « F BLOCK » 

IH HA — ~ Ci HA - - * FIB + C2 HA— 

* LIMITING THE SHUNT CURRENT 

IF < IH HA — . LE. C. 0 > IH HA— - 0.0 

IF ( IH HA — . GT. IMAHA — > IH HA— IMAHA— 

* IH HA — -■-=> I SHUNT 

MACRO RESUME SORT 
END OF MACRO 

tt#**#*##*##*#####*###*##*#**###*######***#*####*-###**##-**!* 

MODEL DESCRIPTION 
LOCATION « 22, HA 
END OF MODEL 
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Physical 
Quantity 

Name I J Description I Units 


VB 

variable 

input bus voltage - 

Volts 

RR1 

parameter 

resistance R1 

Ohms 

RR2 

parameter 

i 

OP-Amp feedback resistance 

Ohms 

CC1 

parameter 

OP-Amp feedback capacitance 

Farads 

RR3 

parameter 

OP-Amp output circuit 
elements 

Ohms 

RR4 

parameter 

(See figure) 

Ohms 

CC2 

parameter 


Farads 

01,02 

parameters 

shunt transistor circuit 
coefficients 


OUTPUT 

Vol 

variable 

output voltage of error 
amplifier 

Volts 

Vcl 

variable * 

See figure 

Volts 

Vo2 

variable 

See figure 

Volts 

Vc2 

variable * 

OP-Amp output circuit 
capacitor voltage 

Volts 

IH 

variable 

shunt current 

Amps 


* state variable 
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MACRO FILE NAME ~ MACROS 

* SHUNT REGULATOR MODEL ( TYPE 1 ) WITH 

* STATE EQUATIONS 
#**#*#*#**##*«****«*##*#*###**#####**##*#******#**###*•»• 

DEFINE MACRO = HB 

MACRO INPUTS =- VC 

RR 1 RR2 CC1 

RR3 RR4 CC2 

Cl C2 

* VB - BUS VOLTAGE 

* RR1 ~ SOURCE RESISTANCE 

* RR2 = OP AMP FEEDBACK RESISTANCE 

* CC1 = OP AMP FEEDBACK CAPACITANCE 

* RR3> RR4> CCS : OP AMP OUTPUT CIRCUIT 

* Cl; C2 : SHUNT TRANSISTOR CIRCUIT COEFFICIENTS 

MACRO OUTPUTS VOl VC1 

V02 VC? IH 

* VOl ~ OUTPUT VOLTAGE OF ERROR AMPLIFIER 

* VC1 » OP AMP FEEDBACK CAPACITOR VOLTAGE 

* V02 = OUTPUT VOLTAGE OF OP AMP OUTPUT CIRCUIT 

* VC2 = OP AMP OUTPUT CIRCUIT CAPACITOR VOLTAGE 

* IH - SHUNT CURRENT 

«***##*##***#***#«£ 

MACRO CODE 
MACRO STOP SORT 

* — — — — - 

* 

* REFERENCE VOLTAGE IS 12. 8 VOLTS 

* SHUNT REF VTG IS 28. 14 

* SENSE VOLTAGE IS (VB SH/2. 1984375) 

it 

it THIS SECTION SIMULATES THE SHUNT ERROR AMPLIFIER 

AA = RR2HB-- / RR1HB — 

BB — i . i- AA 

DD = -AA * (VB HB — /2. 19S4375) + BB * 12.8 + VC1HB — 

it 

* OP AMP VOLTAGE IS LIMITED TO WITHIN 1.5 VOLTS OF 

* SUPPLY VOLTAGES ( VB SH AND GROUND) 

£- 

BB = VB HB — -1. 5 
CC = 1 . 5 

I F ( DD . LT. CC ) DD = CC 
IFU>D . GT. BB ) DD = BB 
VO 1 HB - — - DD 

it 

EE ss VOl HB — - 12. 8 - VC1HB — 

MACRO DERIVATIVE; VC 1HB— = EE / (CC1HB— # RR2HB — ) 

* MODEL OF CIRCUIT AT THE OUTPUT OF CONTROL AMPLIFIER 

<- 

A A = VB HB — - VC2HB — 
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V02HB — - A A 

BB = CC2.HB — * RR3HB — 

FIC2 = AA * < RR3HB — + RR4HB — - ) / ( BB * RR4HB — ) 
S< - VOIHB — / BB - VB HB-- /<CC2HB — *RR4HB — ) 

MACRO DERIVATIVE* VC2HB— ~ FIC2 
«■ 

* MODEL OF SHUNT TRANSISTOR DRIVER CIRCUIT 

K- 

VOC = VOSHB— ~ 

VG = VB HE — - i VOC + 1. 4) 

IF < VG . LT. O. > VG = O. 

FIB = VG / < 66. 5 * 12. ) 

IF ( FIB . LT. 1. 4E-3 > FIB “ 0. 

■* 

* MODEL OF SHUNT ELEMENTS 

* 

IH HB— »■ Cl HE — * FIB + C2 HB— 

IF < IH HB— . IF. 0. 0) IH HB — - 0. 0 
IF < IH HB — . GT. 40. ) IH HB— - 40. 

* IH HB — *««> I SHUNT 

MACRO RESUME SORT 
END OF MACRO 

MODEL DESCRIPTION f ' 

LOCATION ® 22, HB 
END OF MODEL 


/ 
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FULL SHUNT (TYPE 2) [FS] 



INPUT 


Physical 

Quantity 

Name 


Description 

units 

VB 

variable 

Bus voltage 

volts 

VR 

parameter 

Reference voltage 

volts 

R1 

parameter 

1 

ohms 

R2 

parameter 

L Refer to Figure 

ohms 

R3 

parameter 

5 

ohms 

GM 

••'parameter 

Forward transconductance of FET 



OUTPUT 











MACRO FILE NAME * MACROS' 

* FULL SHUNT REGULATOR MODEL (TYPE 2) * 

*-£-*##*##*#****-**#*#**•»■##**###**•#*■#•##*######**•»■ -a-#* 

DEFINE MACRO - FS 

MACRO INPUTS = VB VR VSA R1 R2 R3 

GM VTH 1L.M 

MACRO OUTPUTS - IH VGS 

**#*»«*^**K- •»*•»■!?■•«••»• •«•■«■ ** 

MACRO CODE 
MACRO STOP SORT 

*• — _______ — — — — — — — — ~ _______ — - — 

* VB = BUS VOLTAGE 

* VR « SHUNT REFERENCE VTG 

* VSA “ OPERATIONAL AMPLIFIER SATURATION VOLTAGE 

* GM - TRAiMSCONDUCTANCE OF MOSFET 

* VTH = THRESHOLD VOLTAGE, V(GST), OF POWER MOSFET 

* ILM - HIGH LIMIT OF SHUNT CURRENT 

* IH = SHUNT CURRENT 

* VGS - GATE-SOURSE VOLTAGE OF MOSFET_ 

VD = VB FS VR FS ~ / - 

GAIN = ( R1 FS — •+• R2 FS — ) / ( Ri FS— * R3 FS — > 
VGSFS — = GAIN * VD / GM FS — + -VTHFS — 

* LIMITING THE OP-AMP OUTPUT VOLTAGE 

IF ( VB FS— . LE. VR FS — ) VGSFS — = 0. 

IF ( VGSFS — . GT. VSAFS— ) VGSFS— =■ VSAFS — 

* MOSFET TRANSFER CHARACTERISTIC 

IF ( VGSFS — . L.T. VTHFS— > IH FS — = O. 

IH FS — =■• GM FS— * ( VGSFS VTHFS — ) 

* LIMITING THE SHUNT CURRENT 

IF < IH FS — . GT. ILMFS — ) IH FS — « ILMFS-- 
IF < IH FS— . LE. 0. > IH FS — « 0. 

*• 

MACRO RESUME SORT 
END OF MACRO 

#####*#•&#■!$•*■•*• •***■•«• v*#***#*#*-#*##«# 

MODEL DESCRIPTION 
LOCATION = 22, FS 
END OF MODEL 
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SOLAR ARRAY/PARTIAL SHUNT [AP] 



IL 

IH 



+> VB 


INPUT 


Physical 

Quantity 

Name 

• 

Description 

units 

IL 

variable 

Load current 

amps 

IH 

variable 

Shunt current 

amps 

ILM 

variable 

» * 

Max. solor array output current 

amps 

V1H 

variable 

Max. output voltage of upper array 

volts 

V1L 

variable 

Low limit of output voltage of upper 
array 

volts 

V2H 

variable 

Max. output boltabe of lower array 

volts 


variable 

Low limit of output voltage of lower 
array 

volts 

R 

parameter 

Filter resistance 

ohms 

L 

parameter 

Filter inductacne 

henries 

C 

parameter 

Filter capacitance 

farads 

NP 

parameter 

No. of parallel solar arrays 


NS1 

parameter 

No. of upper serial cells 


NS2 

parameter 

No* of lower serial cells 
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Circuit Model of Solar Array/Partial Shunt [AP] 




ORIGINAL PAQ5 13 
OF POOR QUALITY 


.MACRO FILE NAME MACROS 
DEFINE MACRO AP 

*******************4Hfr*** * ********** ************ 

** SOLAR ARRAY MODEL FOR PARTIAL SHUNT ****** 

ft********************************************** 

*** 


MACRO INPUTS ~ IL TM LL I L.M 

R L. C NP 

NS1 NS2 VI H VI L V2H V2L 

* IL * LOAD CURRENT 

* IH *= SHUNT CURRENT 

* LL » ILLUMINATION LEVEL 

* 1LM « MAXIMUM CURRENT OF SOLAR ARRAYS 

* R» L < C» : CAIILF IMPEDANCE 

* NP = NO. OF PARALLEL ARRAYS 

* MSI a NO. OF UPPER SERIAL CELLS 

* NS2 « NO. OF LOWER SERIAL CELLS 

* V.1H = MAXIMUM OUTPUT VOLTAGE OF UPPER SOLAR ARRAY 

* VI L * LOW LIMIT OF OUTPUT VOLTAGE OF UPPER SOLAR ARRAY 

* V2H a MAXIMUM OUTPUT VOLTAGE OF LOWER SOLAR ARRAY 

* V2L = LOW LIMIT OF OUTPUT VOLTAGE OF LOWER SOLAR ARRAY 

MACRO OUTPUTS = VOi V02 II.. J. 12 

VP 


* VOI a OUTPUT VOLTAGE OF 

* V02 a OUTPUT VOLTAGE OF 

* I Li = UPPER SOLAR ARRAY 

* 12 a LOWER SOLAR ARRAY 

* VB a BUS VOLTAGE 

MACRO CODE 
MACRO STOP SORT 


UPPER SOLAR ARRAY 
LOWER SOLAR ARRAY 
OUTPUT CURRENT 
OUTPUT CURRENT 


*** SOLAR CELL PARAMETERS **«* 

•* 

* RS a INTERNAL SERIES RESISTANCE 

* RSy= INTERN At SHUNT RESISTANCE 

* XIO* REVERSE SATURATION CURRENT 

* TN -- NOMINAL.. TEMPERATURE 

* Q a ELECTRON CHARGE 

* XK a BOLTZMANN CONSTANT 

* XI©* LIGHT-GENERATED CURRENT 

* VOC“ OPEN CIRCUIT VOLTAGE 

K- 

RS a . A 2 
RSH .w 250. 

XIG a .14115 
XIO a 4. je69E-.ll 
A =*:. . 767 
TN * 30.1. 

Q = 1.602E-. -19 
XK = 1. 38J E-23 
VOC = . 5512 
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* XKO ~ Q / < XK * A * TN ) 

XKO-37. H 

■ft 

.* EFFECT OF ILLUMINATION CHANGE 
X I G = X 1 G * LL AP — 

Cl « < 1. F RS / RSH ) 

C21 = NP AP— / < NS 1 AP — * RSH > 

C.22 » NP AP— / < NS2AP — - * RSH ) 

C3 — -NP AP— ■ ft XIG 
All ~ XKO / NSIAP — 

A 1.2 - XKO / NS2AP — 

A2 - XKO «• RS / NP AP — 

IF ( TIME . NE. 0. ) GOTO Ff+5 

ft 

ft INITIAL GUESS FOR SOLAR ARRAY OUTPUT VOLTAGE 

VO LAP — •• < NSIAP — / < NSIAP— + NS2AP — 5 ) 

V02AP- — VP AP— - VO .1 AP — 

•K- ft *8* -JS- ir •& ♦&* ■&•*& *J£ *£• •& ■£* # *8* 

ft ft LIMITING OF SOLAR ARRAY OUTPUT CURRENT BY ft 

ftft MAXIMUM SOLAR ARRAY CURRENT -ft 

ftft I LX * LIMITED SOLAR ARRAY OUTPUT CURRENT ft 

ftftftftftftftftftftftftftftft-ftftftftftftftftft'ftftftftftftftftftftftft-ftftftftftftftftftftftftftft 

fff5 CONTINUE 

Fit X * JL 1 AR- 
IF ( FI LX . GT. ILMAP — > FI LX =•• ILMAP — 

IF ( FIL X . L.T . 0. ) FH...X .« 0. 


* ft * * * ft * ft * * ft ft * * * * * -ft # ft ft # * ft ft * ft ft ft # ft ft ft * ft ft ft ft ft ft ft * ft ft ft ft ft ft ft ft ft 

ftftft NEWTON ITERATION FOR UPPER ARRAYS *********** 

10 CONTINUE 

VOIP - VO JAP — 

FV = Cl ft FI LX + €21 ft V01 AP — F C3 F NP AP — 

& ft XI 0 ft DEXP < All ft VO 1 AP — + A2 * FI LX 

PFV - CP 1, F All ft NP AP — * XIO ft DEXP ( All ft 
+ A 2 ft ! I L. X ' 

V01AP— -- VO 1 AP — - FV / DFV 

2Z1 = < VO 1 AP — - VOIP ) / VO 3 AP — 

IF < DABS < ZZ1 ) , LE. 1 . E-4 ) GOTO +-F-FSO 
GO TO ff MO 


ft LIMIT THE UPPER SOLAR 
F+F20 IF ( V01AP -'.GT.. 

IF ( VO 1 AP — . L.T. 

12 AP — - IL 1AP- — 
IF ( 12 Ar- - . GT. 
IF ( 12 AP — . LT. 


ARRAY OUTPUT VOLTAGE 
V 3 HAP — ) VO 3 AP— =■- V 1 HAP - 
VI LAP— ) VO 3 AP— = VI L AP- 


+ J.H AP — 

ILMAP — ) 12 AP — 
0. ) 12 AP — * 0. 


ILMAP — 


****************** ********************************** 
ftftft NEWTON ITERATION FOR LOWER ARRAYS *********** 
**************************************************** 

+++30 continue: 

VU2P i- V02AP — 


VB AP 


VO 3. AP -- 
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OfttGSMsL 
OF POOR 



FV * Cl * 12 AP — + C2P * VGPAP — . +• CP + WP AP— - 

St * XIO * DEXP ( A 12 * VOPAP’ — + A2 * 12 AP - - ) 

DFV ~ C22 +■ A .1 2 * NP AP— * XIO * DEXP < A12 * V02AP 
& * AP * 12 AP — > 

V02AP-- -• V02AP — - FV / DFV 

222 - ( V02AP— - V02P > / V02AP — 

IF < DABS ( 2Z2 ) . LE. l.E-4 ) GOTO ++ MO 
GO TO h i-30 

* LIMIT THE LOWER SOLAR ARRAY OUTPUT VOLTAGE 
+ *•>40 IF < V02AP — - . GT. V2HAP — - ) V02AP — * V2HAP — 

IF ( V02AP — . L.T. V2LAP— ) V02AP — *= V2LAP — 

VL = VO 1 AP — + V02AP — - IilAP — * R AP — - VB AP — 

MACRO DERIVATIVES, V33 AP — “ '< ILiAP • II AP — ) / C AP 

MACRO DERIVATIVES, I HAP— *■ VL / L AP — 

MACRO RESUME SOFT! 

END OF MACRO 

MODEL DESCR I FT I ON 
LOCATION == 22, AP 
END OF MODEL 
PRINT 
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III. Example of Sysytem- Level Modeling and Simulation 


1. Introduction 

In this section, one example of a spacecraft power system 
is modeled and simulated to show users how to apply the macro 
component models. In addition, the general procedures of 
system-level modeling and simulation are briefly explained 
for the novice. 

The EASY5 commands used in this example are the commands 
used in an interactive mode. 

2. System Model Generation and Simulation 

Step 1. Draw a Block Diagram of the System 

An engineering block diagram of the system to be modeled 
may be drawn to identify the components and their intercon- 
nections of the system. The block diagram of a solar array 
switching system with a buck converter is shown Fig.l. 

Step 2. Write Model Description Statements 

Once required component models and their interconnections 
are identified, appropriate macro components and/or standard 
components are collected from the macro component model li- 
brary and the EASY5 standard component library. 

As shown in Fig. 2, the EASY5 system model [NASA2 .MOD] can 
be described with simple mnemonic statements. 
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Fig. 1 Engineering block diagram of solar array switching system 












The LOCATION command phrase indicates the start of a new 
component in the system model. This command must be followed 
by a numeric value phrase that specifies the location of the 
component on the model schematic. In the example of Fig. 2, 
the location number of the solar array switching unit [AS] 
is 1 and the shunt regulator [FS] is 23, etc. The location 
number phrase is followed by the name of the component at 
that location. A LOCATION statement must be given only once 
for each statement. This means that once a LOCATION state- 
ment is started for a component, the complete description of 
that component must be given. The location of each component 
model should be arranged so that the interconnections among 
the components can be clearly visualized in the EASY5 Model 
Generation program generated schematic diagram. 

Each component model is described with an 'INPUTS' phrase. 
The input component name must be supplemented by the name of 
the particular output quantity that is to provide the input. 
As an example, consider the load component [LO] in Fig. 2. 
Since the output voltage (V2) of the buck-converter power 
stage [BC] is to be the input voltage (VI) to the load com- 
ponent, the following statement indicates to the program that 
the output of [BC], V2, is to be used as the input to the load 
component, [LO]: 

LOCAT I 0N= 10, LO, INPUTS=BC ( V2=V1 ) 
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The input and output connections between the component models 
are made by using the port variables shown with the right and 
left side arrow of the box in the input-output list of each 
component model. 

In model connection an 'implicit loop' should always be 
avoided. The 'implicit loop' is formed whenever the blocks 
are connected with all non-state variables and the connection 
comprises a loop. If a system consists of an implicit loop, 
the Model Generation program will not generate a model, in- 
stead, it will give the error message showing where in the 
system an implicit loop occurs. 

Step 3. Run the EASY5 Model Generation Program 

Whenever a new macro component model is developed, the 
macro model should be compiled and linked by the EASY5 model 
generation program using the command as follows: 

EASY5 AS AS. MOD 


Once the macro cbmponent model is compiled and linked, the 
executable file for the model is stored in the user's perma- 
nent file ' MACROS . DAT ' for subsequent use. After all the 
component models are compiled and linked, the system model 
[NASA2.MOD] should be also compiled and linked by the EASY5 
Model Generation program. With the file ' NASA2 . MOD ' , the 
EASY5 Model Generation program will make all the connections 
between the component blocks according to the Model De- 
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scription data. The command for the EASY5 Model Generation 
program is as follows: 

EASY5 NASA2 NASA2 . MOD 

With the above command, the EASY5 program creates a new file, 

’ NASA2.MGL', which contains a FORTRAN listing of component 
model programs (Fig. 4), a schematic diagram (Fig. 3) and input 
data requirements list (Fig. 5). 

Step 4. Write and Run Analysis Program 

In order to conduct a specific analysis, parameter values, 
initial conditions and any other necessary data in the input 
data requirements list must be supplied. Integration con- 
trols, output data formats and variable names for plot are 
specified as needed. The analysis program of the system is 
given in Fig. 6. 

The command for running the analysis program is 


EASY5 NASA2 * NASA2.ANC 

With this command, the EASY5 creates new files such as 
NASA2. APL, NASA2.PPT and NASA2.RPD. 


Step 5. Analyze the results 
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LIST OF MACROS MACRO COMPONENTS 

COMPONENT NO. 5 NAME = AS 
INPUTS OUTPUTS 


NAME PORT DIM 

NAME 

PORT DIM 

IX 

VO 


IN 

NPP 


DL 

10 


R 

VC1 

STATE 

L 

JL • 

STATE 

C 

R 1 
C 3 
NP 

VB 

STATE 

NS 



MNP 




MACRO CODE 

MACRO STOP SORT 
RS-. 42 
RSH=250. 

X I Q~. 14 i J 5 
X 10=4. 1369E-U 
A—. 967 
TE~ 30 J . 

Q- 1 . 602'E -19 
XK-- 1 . 33 J E— 23 
VOC« 5512 
XK0=39. 3 

IF (TIME. »T. 0. ) GOTO hHO 
NPPAS" — NP AS — 

FT I ME =- 0. 

:-+ t-10 CON'flNUE 

IF ( IH AS — . GE. 5. .AND. NPPAS — . GT. 64 ) GOTO *■ *■ (-20 

IF ( IH AS — - . LT. 2. . AND. NPPAS— . L.T. 324 ) GOTO *■++ 30 

GOTO +W-40 

1-++20 IF < TIME . LT. FTIME-f-DL AS — ) GOTO H-HO 

NPPAS— =• NPPAS MNP AS - 

FT I ME - TIME 
GOTO KM-4Q 

+++30 IF ( TIME . LT. FT I ME H)L AS — ) GOTO ++ MO 
NPPAS— - NPPAS — + NNPAS — 

WRITE <6, «•) IH AS — 

FT I ME =- TIME 
GOTO C 

J !-:-40 CONTINUE 

LL. AS— *- 1. 

X I G = X I G * LL AG — 

Cl = f 1. + RS/RSH ) 

C?~ NPPAS — / < NS AS — * RSH ) 

C3= -NPPAS— *• XJ.G 
AT- XKO / NO AS— - 
A2~ XKO * RS / NPPAS-" — 

IF ( TIME . ME. 0.) GOTO 1- f-50 
VO AS - — r VB AS — 

++*•50 CONTINUE 

VOP=VO AS 


Fig. 4 Part of FORTRAN listing of [AS ] 
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INPUT DATA REQUIREMENTS LIST 
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Fig. 5 Input data requirements list 



OF 


T 1 1 L Jr. -NAG A < LARGE SIGNAL ) 

»/• -v- -Jf* •yr 4r» >t*- ^ •£*• >}> «w» -y“ *$J* ■>>- ^ v* ^ ■>- 45* •fj- "JJ* •is* -<S- •$> *}J* *J5* *j>* *x* *k* •£** ir *x* ^ v- >* •£“ •><■ *i<* ■fj’ ■$<* •>?■ ■Jf* ■x* *ic 

s-TITLE — SOLAR ARRAY SWITCHING SYSTEM WITH 

* SHUNT REGULATOR AND A BUCK CONVERTER LOAD 

■»**«•■»*•»•■«•«••»••«••«•«•«•«•«•*•»*«•#•«••«••«•'»«•«•*•«••»• •S’**#********#****### 

PARAMETER VALUES 

H AS « .001, C AS * 5E-5, L AS = IE -6 

Ci AS — 1. E-3i Ri As — i. O 

NP AS = 324, NS AS = 08, MNP AS -20, DL AS *> .001 

# #• -K- X- * «• -S- * * * £- -K- «■ •£- * * * K- * •«• * Vr fi- •& 

I N I T I AL CONDI I I ONS 

VB AS = 28. 1 55, IL AS = 29.479, VC I AS = 28. 184 
## »«■*#** 

P AR AMETER VALUES 

VR FS « 28. 14, VSAFS = 15. . 

Rl PS - 12400, R2 FS = 2. 2E6, R3 FS « 1. 784 19 

GM FS « 4, VTHFS - 3, ILMFS ~ 45 

C l MC~ 1 , C-2 MC- 1 , C3 MC=1, C4 MC-0 


INITIAL CONDITIONS 

II FI = 27.935, V2 FI - 28. 

PARAMETER VALUES 

R FI=0. 05, L FI-3E-6, C FI= I000L-6 

**JBG (POWER STAGE) 

INITIAL CONDI'! IONS 
IL. BC~9. 47, VC PC- 20 
P AR AMETER VALUES 
L BC- 1 E-4 , C BC~ 4E -4 
RL BC-5E-2, RC BC= IE -1 


**VTG LOOP C0MPENSA1ER 
K ZP=^. 3, ER ZP- 6 
WM ZP--4E3 

WZ1ZP-5E4, WZ2ZP— 1 . 25E2 
WP ZP— 1E6 


K-« PWM 

71 WM=20E-6, VP WM- 6, VQ WM*. 5 
ER WM--&, CICWN-O, SCMWM-0 

** LOAD 

RA LO— 0. 65, RIB LQ-. 6, TC L0-100E-3 
L L.0- 1 E~6, C LO- IE -6 
I N I T I AL CQNDIT I ONS 
XI ZP-34400, X2 2P= . 73, VE ZP*. 3 
IL L0=10, VC L0= 20 

PRINTER PLOTS 
ONLINE PLOTS 
I NT MODE-4 

TMAX— 3E-3, T INC-- PE -7 
PRATED 100, 0UTRA‘!E=100 

SIMULATE 

XIC-X 

P AR AMET ER VALUE S 

RA LO - 0. 65, RB .1.0 = . 76, TC LO = PE -3 
DISPLAY! 

VP BC 

I NT MODE-4 
DISPLAY 5 
VB AS 

DISPLAY2 (OVERPLOT) 

IH FS, 10 AS, I 1 FI 
DISPLAYS, NPPAS 
TMAX---6E--3, T INC'- PE- 7 
PRATE-50, OUTRAIE- 100 
SIMULATE 


Fig. 6 User's Parameter Values and Analysis Commands [NASA2.ANC] 
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The printed data of various variables and parameter values 
at a specified time interval are found in NASA2 . APL ( Fig . 7 ) . 
The plotting data for a line printer is in NASA2 . PPT(Fig. 8) , 
and NASA2.RPD contains a plotting data for a graphics termi- 
nal. The simulation result shown in Fig. 9 was drawn from the 
data in NASA2.RPD. 
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/*/*/*/ SIMULATION ANALYSIS /*/*/*/ 
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Fig. 7 Part of [NASA2.APL] 
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Fig. 9 Simulation results 
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